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^ ■ Abstract: This paper is devoted to the study of various aspects of projectable F{R) 

■ Hofava-Lifshitz (HL) gravity. We show that some versions of F{R) HL gravity may have 

stable de Sitter solution and unstable flat space solution. In this case, the problem of 
scalar graviton does not appear because flat space is not vacuum state. Generalizing the 
C/(l) HL theory proposed in arXiv: 1007.2410, we formulate U{1) extension of scalar theory 
and of F{R) Hofava-Lifshitz gravity. The Hamiltonian approach for such the theory is 
developed in full detail. It is demonstrated that its Hamiltonian structure is the same as 
for the non-relativistic covariant HL gravity. The spectrum analysis performed around flat 
background indicates towards the consistency of the theory because it contains graviton 
with only transverse polarization. Finally, we analyze the spatially-flat FRW equations for 
U{\) invariant F{R) Hofava-Lifshitz gravity. 
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1. Introduction 

In 2009 Petr Hofava formulated new proposal of quantum theory of gravity that is power 
counting renormalizable |||]. This theory is now known as Hofava-Lifshitz gravity (HL 
gravity). It was also expected that this theory reduces to General Relativity in the infrared 
(IR) limit. The HL gravity is a new and intriguing formulation of gravity as a theory with 
reduced amount of symmetries and this fact leads to remarkable new phenomena ^. 

The HL gravity is based on an idea that the Lorentz symmetry is restored in IR limit 
of given theory and can be absent at high energy regime of given theory. Explicitly, Hofava 
considered systems whose scaling at short distances exhibits a strong anisotropy between 
space and time, 

x' = lx, t' = lH. (1.1) 

In (D + 1) dimensional space-time in order to have power counting renormalizable theory 
requires that z > D. It turns out however that the symmetry group of given theory 

^For review and extensive list of references, see lal. 
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is reduced from the full diffeomorphism invariance of General Relativity to the foliation 
preserving diffeomorphism 



x'' = x' + C{t,x) , t' = t + f{t). (1.2) 

The common property of all modified theories of gravity is that whenever the group of 
symmetries is restricted (as for example 1.2) one more degree of freedom appears that is 



a spin-0 graviton. An existence of this mode could be dangerous for all these theories (for 
review, see Q). For example, in order to have the theory compatible with observations one 
has to demand that this scalar mode decouples in the IR regime. Unfortunately, it seems 
that this might not be the case. It was shown that the spin-0 mode is not stable in the 
original version of the HL theory [Q] as well as in SVW generalization Q . Note that in both 
of these two versions, it was all assumed the projectability condition that means that the 
lapse function depends on t only. This assumption has a fundamental consequence for the 
formulation of the theory since there is no local form of the Hamiltonian constraint but the 
only global one. Even if these instabilities indicate to problems with the projectable version 
of the HL theory it turns out that this is not the end of the whole story. Explicitly, these 
instabilities are all found around the Minkowski background. Recently, it was indicated 
that the de Sitter space-time is stable in the SVW setup |^] and hence it seems to be 
reasonable to consider de Sitter background as the natural vacuum of projectable version 
of the HL gravity. This may be especially important for the theories with unstable flat 
space solution. 

On the other hand there is the second version of the HL gravity where the projectability 
condition is not imposed so that N = N{x,t). Properties of such theory were extensively 
studied in 0. It was shown recently in ||8| that so called healthy extended version of such 
theory could really be an interesting candidate for the quantum theory of reality without 
ghosts and without strong coupling problem despite its unusual Hamiltonian structure Q . 
Nevertheless, such theory is not free from its own internal problems. 



Recently Hofava and Melby-Thompson [10| proposed very interesting way to eliminate 



the spin-0 graviton. They considered the projectable version of the HL gravity together 
with extension of the foliation preserving diffeomorphism to include a local U{1) symmetry. 
The resulting theory is then called as non-relativistic covariant theory of gravity ^. It was 



argued there |1C] that the presence of this new symmetry forces the coupling constant A to 



be equal to one. However, this result was questioned in |l^] (see also ||12|, |lj]) where an al 



ternative formulation of non-relativistic general covariant theory of gravity was presented. 



Furthermore, it was shown in |10, 13] that the presence of this new symmetry implies that 
the spin-0 graviton becomes non-propagating and the spectrum of the linear fluctuations 
around the background solution coincides with the fluctuation spectrum of General Rel- 
ativity. This construction was also extended to the case of RFDiff invariant HL gravities 
||8|, |l^ in [14 1 where it was shown that the number of physical degrees of freedom coincides 



with the number of physical degrees of freedom in General Relativity. 



^This theory was also studied in [|l|, ^ |l3[ |l|] . 
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The goal of this paper is to extend above construction to the case of F{R) HL gravities. 
These models were developed in series of papers 

H 3. HL gravity can be considered 



as natural generalization of covariant F{R) gravity. Current interest to F{R) gravity is 
caused by several important reasons. First of all, it is known that such theory may give 
the unified description of the early-time infiation and late-time acceleration (for a review, 
see ]l7| , |2l|].) Moreover, the whole sequence of the universe evolution epochs: Inflation, 
radiation/matter dominance and dark energy may be obtained within such theory. The 
remaining freedom in the choice of F{R) function could be used for fltting the theory 
with observational data. Second, it is known that higher derivatives gravity (like R^- 
gravity, for a review, see ||2^ ) has better ultraviolet behavior than conventional General 
Relativity. Third, modifled gravity is pretending also to be the gravitational alternative 
for Dark Matter. Fourth, it is expected that consistent quantum gravity emerging from 
string/M-theory should be different from General Relativity. Hence, it should be modified 
by fundamental theory. Of course, all these reasons remain to be the same also for the HL 
gravity. Additionally, it is expected that such modification may be helpful for resolution 
of internal inconsistency problems of the HL theory. Indeed, we will present the example 
of F{R) HL gravity which has stable de Sitter solution but unstable fiat space solution. 
In such a case, the original scalar graviton problem formally disappears because fiat space 
is not vacuum state. Hence, there is no sense to study propagators structure around fiat 
space. The complete propagators structure should be investigated around de Sitter solution 
which seems to be the candidate for vacuum space. 

The paper is organized as follows. In the next section, we briefly review the con- 
struction of F{R) HL gravity. Section ^ is devoted to study the de Sitter solutions and 
flat space solutions in such theory. Their stability is analyzed and it is shown that some 
versions of the theory may have stable de Sitter but unstable flat space solution. Clearly, 
this is indication that for such theories the appearance of scalar graviton is not a problem 
due to fact that flat space is not vacuum solution. The whole spectrum analysis should 
be developed around de Sitter vacuum. U{1) extension of F{R) HL gravity as well as of 
scalar HL theory is given in section |^. To check the consistency of such construction, two 
alternative approaches to such extension are proposed. The Hamiltonian structure of U{1) 
invariant F(R) HL gravity is carefully investigated in section ^. It is demonstrated that 
its Hamiltonian structure is the same as in non-relativistic covariant HL gravity. We also 
argue on the general grounds of the Hamiltonian formalism of constrained system that the 
number of physical degrees of freedom coincides with the number of physical degrees of 
freedom in F{R) gravity despite of the absence of the local Hamiltonian constraint. It is 
shown that for the special case F{x) = x the Hamiltonian structure of U{1) invariant F{R) 
gravity coincides with the Hamiltonian structure of non-relativistic covariant HL gravity 
found in which can be considered as a nice check of our analysis. The fluctuations 
around flat background for U{1) invariant F[R) are studied in section ^. It is explicitly 
demonstrated that perturbations spectrum contains the graviton with only transverse po- 
larization. This indicates that scalar graviton problem may be solved within such theory. 



For further study in given direction, see K hg], and for review, see 



-3- 



Section ^ is devoted to demonstration that spatially- flat FRW solutions of U{1) invariant 
F{R) gravity coincide with the ones for the same theory without U{1) invariance. Some 
summary and outlook are given in last section. 



2. Brief Review of F{R) HL Gravity 

In this section we give a brief review of F{R) HL gravity (for more extensive review, 
see This theory is naturally formulated in ADM formulation of gravity Let 

us consider D + 1 dimensional manifold A4 with the coordinates x^, = 0, ■ ■ ■ , D) and 
where = (t, x), x = (x^, • • • , x^). We assume that this space-time is endowed with the 
metric g^u{x^) with signature (— , +,•••,+). Suppose that M. can be foliated by a family of 
space-like surfaces St defined by t = x^. Let gij, {i, j = 1, ■ ■ ■ , D) denotes the metric on St 
with inverse g^^ so that gijg^^ = 5^. We further introduce the operator Vj that is covariant 
derivative defined with the metric gij. We introduce the future-pointing unit normal vector 
to the surface St. In ADM variables one has = \/ —g^^ , = — \/—g^'^- We also 
define the lapse function N = 1/y^—g'^^ and the shift function A^' = —if^/g^^- In terms of 
these variables the components of the metric g^y are written as 

500 = -A^^ + Nig'-Wj , goi = Ni , gij = gij , 

9 Ar2 ' ^ Af2 ' ^ ^ A^2 • 

Then it is easy to see that 

v/-det5 = A^Vdit^. (2.1) 
The extrinsic derivative is defined as 

Kij = ^ (dtgij - ^iNj - V,Ni) . (2.2) 

It is well-known that the components of the Riemann tensor can be written in terms of 
ADM variables. For example, in case of Riemann curvature we have 

{D+i)j^ = K'^K,, -K^ + R+ -^a^ ( x/^n^x) - -^di {^g'^djN) , (2.3) 

where R is D— dimensional curvature. The general formulation of Hofava-Lifshitz F{R) 
gravity was presented in series of papers in [16] ^. This construction is based on the 



modification of the relation (2.3). In fact, the action introduced in [^] takes the form 



Sf{r) = dtd^'x^NFik) , (2.4) 



where 



R = K.jG'^'^'Kki + y^gn^K) - -^d, {Vg9''d,N) - V(ff) , (2.5) 



For further study in given direction, see K hg] , and for review, see 
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where /i is constant, K = Kijg^^ and where the generahzed de Witt metric Q^^^^ is defined 
as 

G'''' = \{g'^g'' + /<7^') - Xg'^g^' , (2.6) 

where A is real constant. Finally V{g) is general function of (^jj and its covariant derivatives. 
We should also note that the special case of F{R) HL gravity where fi = was introduced 
in An important drawback of the case /i = is that it cannot lead to the FRW 

cosmological equations directly. The FRW equations may be obtained only as the limit 
fi = from general theory. In fact, due to the absence of the derivative terms in (^^) this 
theory cannot flow to F{R) gravity action in IR. 



We conclude this section with the remark that the action (2.4) is invariant under 



restricted group of symmetries with is foliation preserving diffeomorphism 

t'-t = fit), x''-x'=eit,x). (2.7) 

Note that under this transformations the metric components transform as 

Nl{x',t') = N,{x,t) - Ni{x,t)f - Nj{x,t)die(.x,t) - g,j{x,t)e{x,t) , 

N'{t') = N{t)-N{t)f, 
g'ij{x',t') = gij{x,t) - gik{x,t)dji^{x,t) - diC^{x,t)gkj{x,t) . 

3. Solutions and their stability in F{R) HL gravity 

In this section, we consider different solutions of standard F{R) Hofava-Lifshitz gravity, 
specially de Sitter and vacuum solutions. The stability of this kind of solutions is studied, 
and it is shown that it depends completely on the choice of function F{R). This suggests 
the way to resolve the scalar graviton problem of Hofava-Lifshitz gravity: Its F{R) version 
may have stable de Sitter vacuum but not flat-space which turns out to be unstable. In 
this situation, the problem is solved simply due to the fact that space flat is not vacuum. 
In order to study the consistency of the theory its spectrum in de Sitter space should be 
investigated what lies beyond the scopes of this work. 

3.1 Stability of de Sitter solutions in F{R) gravity 

Let us consider the stability of the de Sitter solution. As dark energy and even inflation 
may be described (in their simplest form) by the de Sitter space, its stability becomes very 
important topic. Especially in the case of inflation, where a graceful exit is needed to enter 
to radiation/matter dominance, de Sitter space should be unstable. In general, standard 
F{R) gravity contains several de Sitter points, which represent critical points (see [p^]). 
This analysis can be extended to F{R) Hofava-Lifshitz gravity. Let us write the first FRW 
equation in F{R) Hofava-Lifshitz gravity 

= F{R) - 6 f(l - 3A + 3n)H^ + ^i/l F'{R) + QfiHRF"{R) - K^m - , (3.1) 
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For a given F{R), de Sitter solution H{t) = Hq, where Hq being a constant, has to satisfy 
the first equation FRW equation ( ^.ID , 



= F{Ro) - 6^2(1 - 3A + 3/i)F'(i?o) , (3.2) 

where C = and it is assumed the absence of any kind of matter. The scalar R is given in 
this case by, 

^0 = 3(1 - 3A + 6fi)H^ . (3.3) 

Then, the positive roots of equation ( |3.2D are the de Sitter points allowed by a particular 
choice of an F{R) function. Assuming de Sitter solution, one can write F{R) around Rq 
as a series, 

F{R) = Fq + F^{R -Rq) + ^{R- RQf + ^{R- RQf + 0{R^) . (3.4) 

Here, the primes denote derivative with respect to R while the subscript means that it 
is evaluated in Rq. Then, we can perturb the solution writing the Hubble parameter as, 

H{t) = Hq + 5{t) . (3.5) 



Using the function F{R) evaluated around a given de Sitter solution (|3.4D, and the per- 
turbed solution (|3.5|) in the first FRW equation (|3.1D, it yields. 



= ^Fo- 3^2(1 -3A + 3/U) 

-3Fo [((1 - 3A)F^ + QFl^Hl{-l + 3A - 6/i)(-l + 3A - 3/x)) 5{t) 
+6FoV-??o(-l + 3A - 3^) j(t) - l2FQ^l^{t)\ . (3.6) 

Here, we have taken the linear approach on 5 and its derivatives. Note that the first two 
terms in the equation (^) can be dropped because of the equation (|3.2|), which is assumed 
to be satisfied. Then, equation (p. 61) can be written in a more convenient way as, 



6{t)=Q. (3.7) 



m + ^oa-3A+9,).^^^ 

2^ 



1 

+ 



12/i^ 



(3A - 1)-^ - QH^{1 - 3A + 6//)(l - 3A + 3/i) 
-^0 



Then, the perturbations of de Sitter solution will depend completely on the model, that 
is on the derivatives of the F{R) function, as well as on the parameters of the theory 
{\,^Ji)- Note that the instability will be large if the term in front of 5{t) is negative, as 
the perturbations will increase exponentially, while if we have a positive frequency, the 
perturbations will behave as a damped harmonic oscillator. During dark energy epoch, 
as the scalar curvature is very small, the IR limit of the theory can be assumed, where 
General Relativity is recovered, and in such a case we have A = /x ~ 1, and the frequency 

F' 

will depend completely on the value of -g^. In order to avoid large instabilities during the 

F' o 

dark energy phase, the condition > 12i?Q has to be imposed. Nevertheless, during the 
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inflationary epoch, where the scalar curvature is large, the IR limit is not a convenient 
approach, and the perturbations will depend also on the values of (A,/i). Although if we 
assume a very small Fq , the first term in the frequency of the equation ( |3.7D will dominate 



F 

and if A > 1/3, the stability of the solution will depend on the sign of -p^, being stable 
when such a coefficient is positive. 

3.2 On flat space solution in F{R) gravity 

Let us now study fiat space solutions in F{R) HL gravity. In this section we restrict to 
the case of 3 + 1 dimensional space-time. A general metric in the ADM decomposition in 
a 3 + 1 space-time is given by, 

ds^ = -N^dt^ + glf {dx' + N'dt){dx^ + N^dt) , (3.8) 

where i,j = 1,2,3, is the so-called lapse variable, and A^* is the shift 3- vector. For fiat 
space the variables from the metric (|3.8| ) take the values, 

A^ = 1 , Aj = and gij = go6ij , (3.9) 

where go is a constant. Then, the scalar curvature R = 0, and so that our theory has fiat 
space solution, the function F{R) has to satisfy, 

F(0) = 0. (3.10) 

Hence, we assume the condition ( |3.1[1| ) is satisfied, otherwise the theory has no fiat space 
solution. We are interested to study the stability of such solutions for a general F{R), by 
perturbing the metric in vacuum (^^), this yields, 

N=l + 5N{t) and gij = goil + Sgit))6^j . (3.11) 

For simplicity, we restrict the study on the time-dependent perturbations (no spatial ones) 
and on the diagonal terms of the metric. Note that, as we are assuming the projectability 
condition, by performing a transformation of the time coordinate, we can always rewrite 
A^ = 1. Then, the study of the perturbations is focused on the spatial components of the 
metric gij, which can be written in a more convenient way as. 



g,j =m{l + j ^ij ~ 5oe^'^*^'*5y • (3.12) 



By inserting ( |3.12| ) in the first field equation, obtained by the variation of the action on 
A^, it yields at lowest order on (5(t), 

12^^ Fl^5{t)6{t) - 6fi^F^'6^{t) + (-1 + 3X)F[)6^{t) = . (3.13) 

Here the derivatives Fq, Fq are evaluated on i? = 0. Then, the perturbation 5 will depend 
completely on the kind of theory assumed. We can study some general cases by imposing 
conditions on the derivatives of F{R). 
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• For the case Fq = 0, it gives 5{t) = 0, such that at lowest order the flat space is 
completely stable for this case. 

• For Fq = and Fq ^ 0, the differential equation ( 3.13| ) has the solution, 

m = Cit[l + ^^+C2 (3.14) 

where Ci^2 are integration constants. Then, for this case, the perturbations grow as 
the power of the time coordinate, and flat space becomes unstable. 

Hence, depending on the theory, the flat space solution will be stable or unstable, 
which becomes very important as it could be used to distinguish the theories or analyze 
their consistency. 

3.3 A simple example 

Let us now discuss a simple example. We consider the function, 

F{R) = kqR + mir , (3.15) 

where ko,i are coupling constants and n > 1. Note that this family of theories satisfies 
the condition ( 3.1C| ). The values of first and second derivatives evaluated in the solution 
depend on the value of n in ( |3.15| ), 

F'(0) = KQ , F"(0) = Kin(n - l)E!^~^ . (3.16) 

Then, we can distinguish between the cases, 

• For n 7^ 2, we have -F"(0) = 0, and by the analysis performed above, it follows that 
flat space is stable 

• For n = 2, we have -F"(0) = ki, and flat space is unstable. 

Hence, we have shown that the stability of solution depends completely on the details of 
the theory. 

We can now analyze the de Sitter solution. Using the equation ( |3.2| ), we can find the 
de Sitter points allowed by the class of theories given in Eq. ( p. 151 ) , 

3„2,o, , (3gg(l - 3A + 6^))" (1 - 3A + 6^ - 2n(l - 3A + 3/.)) 
-Fo(3A-l).o + 2(1 - 3A + 6^) = °- ^^'^'^ 



Resolving the Eq. (3.17), de Sitter solutions are obtained. For simplicity, let us consider 
n = 2, in such a case the equation ( p. 17 ) has two roots for Hq given by. 



3Vki(1 - 3A + 6/i)(-l + 3A - 2ii) 

As we are interested in de Sitter points, we just consider the positive root in ( p. 18 ). Then, 
the stability of such de Sitter point can be analyzed by studying the derivatives of the 
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pi 

function F{R) evaluated in Hq. The stability will depend on the value of j^, which for 







this case yields, 

^0 _ '^0 _^^2Hl. (3.19) 



In the IR limit of the theory (A — )• 1, /x — )• 1), the condition for the stability of de Sitter 
points > IS-fTg is clearly satisfied by ( 3.191 ). Even in the non IR limit, Fq = 2ki and 



- 

assuming ki <^ 1, we have that, 

^ = 3Hi{l-3X + 6^^) + ^. (3.20) 

As this term is positive, we have that the instabilities will oscillate and be damped, such 
that the de Sitter point becomes stable. 

Thus, we presented the example of the F{R) theory where flat space is unstable solution 
and de Sitter space is stable solution. The problem of scalar graviton does not appear in 
this theory because one has to analyze the spectrum of theory around de Sitter space 
which is real vacuum. Indeed, flat space is not stable and cannot be considered as the 
vacuum solution. Of course, deeper analysis of de Sitter spectrum structure of the theory 
is necessary. Nevertheless, as we see already standard F{R) gravity suggests the way to 
resolve the pathologies which are well-known in Hofava-Lifshitz gravity. 

4. U{1) Invariant F(R) Hofava-Lifshitz Gravity 

Our goal is to see whether it is possible to extend the gauge symmetries for above action 



as in [10 1 . As the first step we introduce two non-dynamical fields A,B and rewrite the 



action (|2.4| ) into the form 

^FiR) = dtd''x^N{B{R -A) + F{A)) . (4.1) 

It is easy to see that solving the equation of motion with respect to A, B this action reduces 
into (p.4|). On the other hand when we perform integration by parts we obtain the action 
in the form 

Sf{R) = dtd^'x (y^NB{Kijg'^^'Kki - V{g) - A) 

+^NF{A) - 2fi^NVnBK + 2fidiB^g'^ djN) , (4.2) 
where we ignored the boundary terms and where 

VnB = ^{dtB-N'diB). (4.3) 
Let us now introduce U{1) symmetry where the shift function transforms as 

5aNi{x, t) = N{x, t)\/ia{x, t) . (4.4) 

It is important to stress that as opposite to the case of pure Hofava-Lifshitz gravity the 
kinetic term is multiplied with B that is space-time dependent and hence it is not possible 
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to perform similar analysis as in |10|. This procedure frequently uses the integration by 
parts and the fact that covariant derivative annihilates metric tensor together with the 
crucial assumption that depends on time only. Now due to the presence of B field we 
have to proceed step by step with the construction of the action invariant under ([4.4|). 
As the first step note that under the kinetic term 5*^™ = I dtd^ x^KijQ^^^'- K^i 
transforms as 

In order to compensate this variation of the action we introduce new scalar field u that 
under (4.4) transforms as 

6ai^{t,x) = a{t,x) (4.5) 
and add to the action following term 

Sl^^ = ^ J dtd''x^NBKijg'^'''V,VjU. (4.6) 

Note that under (4.5) this term transforms as 

6aSi^'> = ^ J dtd'^x^NBK.jG'^^'Vkyia-^ j dtd^'x^NBViVjag'^'^'VkViu. 

so that we add the second term into the action 

5(2) = ^ j dtd^x^NBV.Vjvg'^^^VfSJjV. (4.7) 

As a result, we find that S^'"" + 5^^^ + 5^,^^ is invariant under (O) and (O). 



As the next step we analyze the variation of the S-kinetic part of the action S^^^'^ 
/ dtd^Xy/gNVnBK under the variation 



,^^5^kin = ^ y dtd^x^aV'iViBK) + ^ y dtd^x^NaViVjig'^VnB) . (4.8) 

We see that in order to cancel this variation it is appropriate to add following expression 
into the action 

S"-^ = -^fi j dtd^x^NuV\V,BK) - ^iJL j dtd^x^Ni^V^Vjlg'^V^B] 
+ j dtd^x^NV^uVkBViV'u . 



Then it is easy to see that S'-^'^™ + 5''^ is invariant under ( [4.4[ ) and ( |4.5| ). Collecting all 
these results we find following F{R) HL action that is invariant under ([4. 4]) and ([4.5|) 



Sf{r) = d^d""^ [^NB{K,,g'^'''Kki - V{g) - A) 

+^NF{A) - 2f,^NVnBK + 2f,d.,B^g'^ djN) 

-2fi j d^xdty/gvV\V,BK)-2iJL j d^ xdt^NvV -sj j{g'^V ^B) 

+2n j dtd^x^NV'uViBV^VjV 

+2 j d^x^NBK^jg'^^^VkViu + J d^x^BV^Vjug'^'^^VkViu. (4. 
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Note that we can write this action in suggestive form 
Sf(R) = dtd^'x (^^NB{{Kij + V ,V ju)G'^''\Kki + Vk^i^) - V[g) - A) 

(4.10) 

or in even more suggestive form by introducing 

Ni = N,- NV,u , K,j = ^ {dtQij - ViNj - VjNi) (4.11) 

so that 

^FiR) = ^J dtd^'x {^^NB{K,^G'i^'Kki - V{g) - A) 

+^NF{A) - 2^^NV^Bg'^K,i + IfxdiB^g'^ B^n) (4.12) 



is formally the same as the F{R) Hofava-Lifshitz gravity action. Clearly this action is 
invariant under arbitrary a = a{t,x). Moreover, such an introduction of U{1) symmetry 
is trivial and does not modify the physical properties of the theory. This is nicely seen 
from the fact that ly appears in the action in the combination with Ni through iVj where v 
plays the role of the Stiickelberg field. In order to get physical content of given symmetry 



we follow [10 1 and p^ ] and introduce following term into action 

'5"'' = ^ / dtd''x^Bg{gij){A - a) , (4.13) 

where 

N ■ 

a = v-N'ViU + —V'ViU, (4.14) 
where X = Note a transforms under a variation as 

a'{t, x) = a{t, x) + a{t, x) — N^{t, x)Via{t, x) . 
Then it is natural to suppose that A transforms under a-variation as 

A'{t, x) = A{t, x) + a{t, x) - N\t, x)Via{t, x) (4.15) 

so that we immediately see that ^ — a is invariant under a-variation. The function Q 
can generally depend on arbitrary combinations of metric g and matter field and we only 



demand that it should be invariant under foliation preserving diffeomorphism (2/7) and 



under (|4.4|) and ([4.5|) . For our purposes it is, however, sufficient to restrict ourselves to 
the models where Q depends on the spatial curvature R only. Now one observes that the 
equation of motion for A implies the constraint 

Bg{R) = (4.16) 
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that for non-zero B implies the condition G{R) = 0. Note that this condition is crucial for 
elimination of the scalar graviton when we study fluctuations around flat background. We 
demonstrate this important result in the next section. 

Finally note that it is possible to integrate out B and A fields from the actions ( 4.12| ) 
and ( 4.13| ) that leads to 



= ^ j dtd^x^NF{R) 



R = K^.g^^'^'Kki - V{g) + {dt - [VgN'R)] + ^g{A - a) . 

(4.17) 

This finishes the construction of U{\) invariant F{R) HL theory action. 
4.1 Lagrangian for the Scalar Field 

Now we extend above analysis to the action for the matter field with the following general 
form of the scalar field action 



'matt 



j dtd^x^NX , (4.18) 



where 

X = -(V„</.)2 + F{g'^di(j}dj<P) . (4.19) 
where F{x) = X + Yln=2 where we defined 

Vn(^=^{dt(l>-N'dicl>). (4.20) 

Note that this general form of the scalar field action is consistent with the anisotropy of 
target space-time as was shown in [24|. 



Now we try to extend above action in order to make it invariant under ( [4.4[ ). Note 
that under this variation the scalar field action (4.18) transforms as 



SaSra^tt = -2 j dtd^X^NV'aVi<i)VnCt> (4.21) 

using 

5aX = 2V'aVi(j)Vn(p . (4.22) 



We compensate the variation ( 4.21| ) by introducing additional term into action 



Sme.tt-u = -2 y dtd^x^vNV\Vi(j)Vri<l)) + j dtd^x^V'vV^uV,(j)Vj(l) (4.23) 
Then the action ( [4.23 ) transforms under ( |4.4| ) and (^]^) as 

6aSrr.^tt-v = 2 j dtd^x^NV'aVi(j)Vn<j) (4.24) 



that compensates the variation ( 4.21| ). 
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In the same way one can analyze more general form of the scalar action 

5matt = -J dtd^x^NKiX) . (4.25) 



In order to find the generalization of the action ( 4.25| ) which is invariant under (4.4) we 



introduce two auxiliary fields C, D and write the action ( ^.25| ) as 

5matt = - J dtd''x^N[K{C) + D{X - C)] . (4.26) 
Clearly this action transforms under (4^) as 

(^a^matt = - j dtd^x^NDd^X = -2 j dtd^ x^N DV' aV iC^V ncj) , (4.27) 
where relation ( |4.22| ) is used. It is easy to see that the variation of the following term 

5'matt-z. = 2 j dtd^x^NDV'vVi<t)Vn(ly + j dtd^x^DV'uV^vVi(INj(i) (4.28) 



compensates the variation ( 4.27 ). Finally note that ( 4.26 ) together with ( [4.28| ) can be 



written in more elegant form 

^matt = - j dtd^x^N[K{C) + D{X -C)] = - j dtd^x^NK{X) , (4.29) 

where 

X = -(v„(/.)2 + F(5^^a,(/>5,-0), 

Vn(t> = ^{dt<t) - N'Vict>) = ^(5t0 - X'Vi + NV'uVicfy . 

In this section we constructed C/(l)-invariant scalar field action in the form which closely 
follows the original construction presented in In the next section more elegant ap- 

proach to the construction of U{1) invariant F{R) HL gravity and the scalar field action 
is given. 

4.2 Alternative Definition of U{1) Invariant F[R) HL gravity 

To begin with we note that under the (local) C/s(l) symmetry, the shift function Ni{x,t) 
and u{x,t) are transformed as 

Ni{x,t) ^ Ni{x,t) + N{t)Via{x,t) , v{x,t) ^ v{x,t) + a{x,t) . (4.30) 

Therefore the combination 

Ni{x, t) = N,{x, t) - N{t)Viiy{x, t) , (4.31) 

is invariant under the transformation of the local Uj]{l) symmetry. Then if we replace Ni 
with Ni, one can always obtain the model with the local U-£{1) symmetry. 
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For example, the extrinsic curvature Kij could be replaced by 

Kij ^ Kij = ^ {dtgij - ViNj - VjNi) = Kij + ^ {ViVj + V^ V,) u . (4.32) 
Then it follows 



2 

In the same way one can deal with following term 

(n^V.n'^ - n^V.nn = {dt {^K) - d., {^N'K)] 

^^{dt {V~9K) - iVgN^K) } . (4.34) 

The above analysis tells that if we define a "curvature" by 

R = K,,K'i - + [dt iVgK) - iVgN'K) } - Vig,,) , (4.35) 



the action 



S = ^ j dtd^x^NF{R) , (4.36) 
2k J 



is invariant under the local Ut,{1) transformation (|4.30 ). Finally, in order to obtain non- 
trivial symmetry we have to add to R the expression g(gij)(A — a). Then the action derived 
here coincides with the action ( [4.17 ). 

Note that the equation obtained by the variation of gives a constraint which kills 
the extra and problematic scalar mode appearing in the original Hofava gravity. For the 
action ( 4.36| ), the equation has the following form: 

= y (ViV, + V,V,) {K'^F' {R)) - XNV^ {KF' {R}) 

-2^V2 {dtF' {R) - N'diF' {R)) 

I ( d , ,. N. 



+^ [j^{V~9F'{R)g) - V^{N^^F'{r)g) - -V^V^{^F' {R)g) j =0. 

(4.37) 

n the same way we can proceed with the action for scalar field (p ( 4.25| ) . Instead of using step 
by step procedure performed in previous section one immediately makes the replacement 
^ = - iVVV that leads to the action (|429|) . 
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5. Hamiltonian Formalism of ^7(1) Invariant F(R) HL Gravity 
Let us again consider an action 

SpiR) = ^ / (ygNB{Ki^g'i'''Kki - V{g) - A) 

+^ j dtd^x^Bg{gij){A-a). (5.1) 

Now we perform the Hamiltonian analysis of given action. Note that the Hamiltonian 
analysis of F{R) HL gravity was performed previously in ||l^, ^ Q and we generalize 
these works to the case of U{1) invariant F{R) HL gravity. 
From ( ^.1| ) we find the conjugate momenta 

= ^Bg'^'''Kkl-^l^VnBg'^ , p'-o, 
K^PB = -2ii^K , PA~0, PA^O, p^ = - — ^Ng. 

From these relations the primary constraints are found 

'^1 : Pa{x) ^ , <^>2:p^{x) + ^^Ng{x), pa{x) ^ , p'{x)^0, pN - . (5.2) 

The Hamiltonian density is obtained in the form 

n = NUt + N'Ui , 

where 

2 2 2 

^B ^BD ^fiD 

- - ^2^9B{V{g) -A)- ^,^gF{A) + ^^d^Vag'^d^B] 



-2uViVjTT'^ + uV'ViB + —g{R)BV'ViU, 



1 

r 

1 



-Hi = -2gaVk^^^ + PB^iB - ■^^Bg{R)ViV , (5.3) 



where 



gijkl = ij{9ik9jl + 9il9jk) - _ ^ 9ij9kl ■ (5.4) 



According to the general analysis of constraint systems ^ we should consider following 
Hamiltonian 

H = j d^xiNTir + N'Tii + v^<^i + + v^pa + Wip' + wnPn) j d^x^g{R)A , 

^For review, see 123]. 
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where f^'^, w^, Wi, and wn are Lagrange multipliers related to corresponding primary 
constraints. Now the requirement of the conservation of the primary constraints pat ~ 0, 
Pi{x) « 0, $1^2(2;) ~ and pa{x) implies following secondary ones 

dt^i = {<^i,H} = -\b^ g = -^{' « , 



dtPN = {PN: ^} = - / ^^^^r = -T « , 

dtPi = {pi,H} = -ni^o, 



= -2Rij7r'^ + -Rn - -BVkVH^) + —(RPB - (1 - D)V,V'pb) ■ (5.7) 



dtPAix) = -^B + ^F'{A) = ^Ga ~ . (5.5) 
Finally we determine the time evolution of <I>2. Using 

{R{x),'n'^{y)] = -R'^{x)5{x - y) + V'V^6{x - y) - g'^VkV'dix - y) 

we find 

dt^2 = {<^2,H}^N^^' , (5.6) 
where the secondary constraint takes the form 

Note that in the calculation of ( p.6| ) we used the fact that 

{p,, H} = -NV'Ui + \y^{^N'g) + Av*V,(V5a) ~ « (5.8) 

using that the right hand of this equation is proportional to <1>(^ or its covariant derivatives 
and hence it vanishes on constraint surface w 0. In fact, we previously simplified the 
equation ( |5.7D where all terms proportional to the constraint were ignored. 

Now one can summarize the results. We have the primary constraints pN ^ 0, Pi ^ 0, 
PA ^0, Pu^ 0, w 0, w 0, Ri 0, T « 0. Then the total Hamiltonian takes the 
form 

H = NT + v^pN 

+ j d^x{N'ni + v'pi + v^pA + v^p^ + v}j<^^^ + vjj<^i^ + v^pA + v^Ga) , 

where we included an expression ^{^A into the expression v\j^{^ . 

Let us analyze the consistency of the secondary constraints with the time development 
of the system for the case F"{A) 7^ 0. To begin with note that one can write Hi as 
Tii = — 2(7jfcV;7r'^' +pBdiB + pydiV — (^[^diV so that it is natural to consider as the smeared 
form of diffeomorphism constraint the following expression 

Ts{M') = j d^xM'{-2gikVi7r''^ + psd.B + p.diu) (5.9) 

that is the generator of spatial diffeomorphism. Then using the fact that the total Hamil- 
tonian is manifestly invariant under spatial diffeomorphism one finds that Tg is preserved 
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on the constraint surface. In case of the constraint p,y, using (^^) and also the fact that 
{p^,^{^} = {pi/,$2^} = it follows that this constraint is preserved during the time 
evolution of the system and also that it is the first class constraint. 

Let us now consider the constraint <I>{^. Its time development is governed by equation 

dM' = {^i',HT}^N^^' + J d^xvh{x){<^i',^'ix)] . (5.10) 

To proceed further we note that from the structure of the constraint <I>(^ and we clearly 
have 

{$f(a;),$^^(y)}/0. (5.11) 

Say differently, together with are the second class constraints. Then the require- 
ment that the right side of the equation ( ^.10| ) has to vanish implies that the Lagrange 
multiplier vjj has to be equal to zero. On the other hand the time evolution of is 
governed by equation 

dtH' = Nm',T}+ I d^x{v^{x){<^2,G{x)}+v}j{x)m',<^{'{x)}) . (5.12) 

Note also that {$2^(a;), G^(y)| / as follows from ( ^!5|) and from ( ^?7|) . On the other 
hand the time evolution of Ga is governed by equation 

dtGA = {GA,HT} = N{GA,T} + I d''x{vh{x){GA,H'{x)}+v^{x){GA,PA{x)}) 

(5.13) 

Since {pAix), GA{y)} = —^/gF"{A)6{x — y) and since v'jj = we see that requirement of 
the vanishing of the right side of ( |5.13| ) determines as a function of canonical variables. 
On the other hand the requirement of the preservation of the constraint pA 

dtPA = {PA, Ht} « J d'^xv^ {pA, Ga{x)} = (5.14) 

implies v'^ = 0. However using this equation in ( |5.12| ) we can again determine the Lagrange 
multiplier vjj as function of canonical variables. At this place we see that the requirement 
of the preservation of the constraints ^•('^2' PA^ Ga determines the corresponding Lagrange 
multipliers v^f, v^, and consequently no new constraints have to be imposed. 

Now we are ready to determine the number of physical degrees of freedom. To do this 
note that there are D{D + 1) gravity phase space variables gij, tt*-', 2D variables Ni,p^, 2 
variables A, pa, 2 variables B, pB, 2 variables A, pA and 2 variables I'jPu- In summary the 
total number of degrees of freedom is A'^o.o.f = + 3-D + 8. On the other hand we have 
D first class constraints Tii ^ 0, D first class constraints pj w 0, 2 first class constraints 
Pu ^ 0, PA ^ ^ and four second class constraints $(^, ${^^, PA, Ga- Then there are 
-^f.c.c. = 21) + 2 first class constraints and A's.c.c. = 4 second class constraints. The number 



of physical degrees of freedom is [25| 



iVD.o.f. - 2iVf.e.c. - A^s.c.c. = {D^-D-2) + l (5.15) 



-17- 



that exactly corresponds to the number of the phase space physical degrees of freedom of 
D + 1 dimensional F{R) gravity. For example, for D = 3 the equation ( 5.15| ) gives 4 phase 



space degrees of freedom corresponding massless graviton and 2 phase degrees of freedom 
corresponding to the scalar. Note also that there is still global Hamiltonian constraint that 
has to be imposed and also that all second class constraints have to be solved. Solving of 
PA = 0, Ga = one can express ^ as a function of at least in principle. Unfortunately 
solving the second class constraints $(^, is very difficult in the full generality. On 
the other hand, it is easy to see that in linearized approximation these constraints can be 
solved as /i = 0, vr = where h is the trace part of the metric fluctuation and vr is its 
conjugate momenta. 

The previous analysis was valid for the case when F"(A) 7^ 0. Let us now discuss the 
second case when F'{A) = 1 ®. For F'{A) = 1 we see that the Poisson bracket between pA 
and Ga is zero. Then the equation ( ^.13 ) implies additional constraint (Note that vjj = 0) 

dtGA = {Ga, Ht] = N {Ga, T] = -^^G'J , (5.16) 



where 



G'i = 7:-{\D-l)^. (5.17) 

Due to the fact that {pa, Ga} = one finds that ^ is the first class constraint. On 
the other hand we have following non-zero Poisson bracket 

{GA{x),G'i{y)] = ^-hR^B{x)5{x-y) (5.18) 

that implies that Ga and G^^ are the second class constraints. The analysis of the remaining 
constraints is the same as in case F"{A) 7^ so that we have following set of constraints. 
There are D first class constraints Hi ^ 0, D first class constraints w 0, 3 first class 
constraints pi^ ^ p_A ^ 0, PA ^ and four second class constraints ${^, ^{^^ , Ga, G^^. 
Then we have A^f.c.c. = 2D + 3 first class constraints and N^.c.c. = 4 second class constraints. 
The first class constraint pA = ^ can be eliminated with the gauge fixing condition A = \. 
Solving the constraint Ga = we obtain B = 1 while solving the constraint = we 
find 

(5.19, 



Inserting this result into T-Lt given in ( ^.3| ) one gets that it takes the form 

that corresponds to the Hamiltonian constraint of non-relativistic covariant HL gravity 
whose explicit form can be found in |14]. In the same way we insert (5.19) into ( |5.7| ) and 



we find 

= -2/?i,-7r*^' + 2—^^ — Rn + VjVV 
2 ^ D\-l XD-1 



® Generally we could have F'{A) — C for constant C. However we show below that in this case B = C 
and consequently C can be eliminated by redefinition of 
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which again coincides with the constraint found in In other words in case when 

F"(A) = the Hamiltonian structure of U{1) invariant F{R) gravity coincides with the 
Hamiltonian structure of non-relativistic covariant HL gravity. This result can be consid- 
ered as a nice confirmation of our procedure. 



6. Study of Fluctuations around Flat Background in U{1) Invariant F{R) 
HL Gravity 

Let us analyze the spectrum of fluctuations in case of U{1) invariant F{R) HL gravity for 
the special case /i = 0. For simplicity we assume that F{R) Hofava-Lifshitz gravity has 
flat space-time as its solution with the background values of the fields 

g^^=6,,, iV(0) = i, nP=0, ^W=0, z.(0)=0. (6.1) 

Note that for the flat background the equation of motion for B and A takes the form 

V(5(0))-^(0) =0, 5W-F'(^W) = (6.2) 

that implies that A^^\ B^^^ are constants. In order to find the spectrum of fluctuations we 
expand all fields up to linear order around this background 

9ij = ^ij + l^^ij > = '^"-j ; N = 1 + KU , 

A = A(°)+Ka, B = B^^^+Kb, A = A^^^ + kA, v = v^^^+Ki>. 

Since n depends on t only, its equation of motion gives one integral constraint. This 
constraint does not affect the number of local degrees of freedom. For that reason it is 
natural to consider the equation of motion for hij, rii and v only. We further decompose 
the field hij and rii into their irreducible components 

hij = Sij + diWj + djWi + ^didj - j^^ijd^^ ^ + j^^ij^^ (6-3) 

where the scalar h = ha is the trace part of hij while Sij is symmetric, traceless and 
transverse 

d'sij = 0,d' = d'^dj (6.4) 

and Wi is transverse 

^'w^ = 0. (6.5) 

In the same way we decompose rii 

Ui = Ui + diC (6.6) 

with Ui transverse d'^Ui = 0. In what follows we fix the spatial diffeomorphism symmetry 
by fixing the gauge 

u;. = , M = . (6.7) 



-19- 



We begin with the equation of motion for Ni 

2VjiBg'^''^KM) + BG{R)ViV = . (6.8) 

and for u 

iv.v.iBNG'^'^'Kki) + (j^iVaBG) - Vii^BN'g) - ^V^V'i^BNg)^ =0. 

Combining these equations and using the equation of motion for A one obtains 

g{R) = (6.9) 

Then 

Bj^g - BN'v^g - ^{BViV'g + 2v,sv'g) = o . (e.io) 



Further, in the linearized approximation the equation of motion ( |6.8| ) takes the form 

- 2Boil - \)d'{dkd^C) + BodkO'^u' + - DX)d'h + 2(1 - X)Bod\djd^D) = (6.11) 

using the fact that g{Ro) = and also d^Sij = d^Sij = and 5"^^ Sji = 0. Let us now focus 
on the solution of the constraint g{R) = in the linearized approximation. Let Rq^^ is 
the solution of the equation of motion and let us consider the perturbation around this 
equation. These perturbations have to obey the equation 

g{Ro + 6R) = g{Ro) + ^dR = %{Ro)5R = . 

aR aR 

We see that in order to eliminate the scalar graviton we have to demand that ^(Ro) 7^ 0. 
To proceed further note that 

5R,^ = i[vf + Vf V'=(°)/i., - V'^W/.., - Vf^vfh] (6.12) 



where V^°^ is the covariant derivative calculated using the background metric glj\ Then 
in the flat background it follows 



5R = -(1 - D)d''dkh 

so that the condition 6R = implies d^d^h = 0. Then h = hit). However, in this case 
the fluctuation mode does not obey the boundary conditions that we implicitly assumed. 
Explicitly we demand that all fluctuations vanish at spatial infinity. For that reason one 
should demand that h = 0. Then it is easy to see that the equation ( |6.10D is trivially 
solved. In the linearized approximation the equation of motion for gij takes the form 

^{sij + ^Sij{l - XD)h - diiij - 2d,djC + 2d'd^D - 2X5'^dkd''i) + 2X6^j^k^''C) 

-— + ^[d'&'i^ - «) + S'^dkd'^iA - a)] = (6.13) 
gij an 
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where 

a = | + ^4a'.. (6.14) 

Note that we have not fixed the U{1) gauge symmetry yet. It turns out that it is natural 
to fix it as 

iy = 0. (6.15) 
Then the trace of the equation ( |6.13|) is equal to 

(1 - XD)dkd''C + 6'^^ - ^(1 - D)^^^'A = . (6.16) 
ogij dR 

Let us again consider the equation of motion ( |6.11| ) and take its 5*. Then using the fact 
that d'^Ui = one gets the condition C = f{t) that again with suitable boundary conditions 
implies C = 0. However then inserting this result in ( |6.11| ) we find 

dkd'^u, = (6.17) 

that also implies Ui = 0. 

To proceed further we now assume that V2 = —R so that = —^dkd^ Sij — ^f^{didj — 
5ijd^)h. Clearly the trace of this equation is proportional to h and hence it implies following 
equation for a 

dkd^A = ^. (6.18) 

Imposing again the requirement that A vanishes at spatial infinity we find that the only 
solution of given equation is ^ = 0. Finally the equation of motion for gij gives following 
result 

Sij+dkd^Sij = Q. (6.19) 

In other words, it is demonstrated that under assumption that U{1) invariant F{R) HL 
gravity has flat space-time as its solution it follows that the perturbative spectrum contains 
the transverse polarization of the graviton only. Clearly, this result may be generalized for 
general version of theory with arbitrary parameter ^u. 

Finally we consider the linearized equations of motion for A and B. In case of A one 

gets 

-b + F"{AQ)a = Q (6.20) 

while in case of B we obtain 

dV 

-^{Ro)5R-a = (6.21) 
Using the fact that SR ~ /i = we get from (|6.21| ) and from ( |6.20| ) 

a = 6 = 0. (6.22) 

In other words there are no fluctuations corresponding to the scalar fields A and B. One 
can compare this situation with the conventional F{R) gravity where the mathematical 
equivalence of the theory with the Brans-Dicke theory implies the existence of propagating 
scalar degrees of freedom. In our case, however, the fact that U{1) invariant F(R) HL 
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gravity is invariant under the foliation preserving difFeomorphism allows us to consider 
theory without kinetic term for B (// = 0). 

Thus, it seems U(l) extension of F(R) HL gravity may lead to solution of the problem 
of scalar graviton. 

7. Cosmological Solutions of U{1) Invariant F{R) HL gravity 



Let us investigate the FRW cosmological solutions for the theory described by action ( 4.36| ). 
Spatially-flat FRW metric is now assumed 

3 

ds'^ = -N^dt^ + a\t)^{dx')\ (7.1) 

i=l 

and we choose the gauge fixing condition for the local ?7s(l) symmetry as follows, 

1^ = 0. (7.2) 



In the FRW metric ( [7.1D , one gets iVj = and R and Kij only depend on the cosmological 
time t. Therefore, the constraint equation ( [4.37| ) can be satisfied trivially. For the metric 
( |7.1| ) and the gauge fixing condition, the scalar R is given by 

- 3(1 - 3A + 6/i)g^ ^ 6fid fH\ 



The second FRW equation can be obtained by varying the action (^[36) with respect to the 
spatial metric gij, which yields 



= F{R) - 2(1 - 3A + 3^) i^H + 3H^j F'{R) - 2(1 - 3X)RF"{R) + 

+2/i f (R) + RF"{R)) + K^pm ■ (7.4) 



Here the contribution from the matter is included. Pm expresses the pressure of a perfect 
fluid that fills the Universe, and = 1. Note that this equation becomes the usual second 
FRW equation for conventional F{R) gravity by setting the constants A = /i = 1. The 
variation with respect to N yields the following global constraint 







j d^x F{R) - 6(1 - 3A + 3yL)H'^ - QfiH + Q^HRF"{R) - K^p^ . (7.5) 



By assuming the ordinary conservation equation for the matter fluid pm+3H{pm+Pm) = 0, 
and integrating Eq. ( |7.4D , 

= F{R) - 6 f(l - 3A + 3p)H^ + pll] F'{R) + 6pHRF"{R) - K2p^ - ^ , (7.6) 

where C is an integration constant, taken to be zero, according to the constraint equation 
( |7.5| ). Eq. ( |7.6| ) corresponds to the first FRW equation. Hence, starting from a given F{R) 
function, and solving Eqs. ( |7.4| ) and (^), FRW cosmological solution can be obtained. 
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Note that the obtained equations (|7.4| ) and ( |7.5D are identical with the corresponding 
equations in the F(i2)-gravity (for recent review of observational aspects of such theory 
see pO| ]) based on the original Hoirava gravity. Hence, U{1) extension does not influence 
the FRW cosmological dynamics. 

Let us consider the theory which admits a de Sitter universe solution. We now neglect 
the matter contribution by putting pm = Pm = 0- Then by assuming H = Hq, Eq. ( [7^ ) 
gives 

= F (3 (1 - 3A + 6/x) H^) - 6 (1 - 3A + 3/x) H^F' (3 (1 - 3A + 6/i) H^) , (7.7) 
as long as the integration constant vanishes (C = 0). We now consider the following model: 

F{R) ocR + pR^ + 7i?3 . (7.8) 

Then Eq. ( |7.7| ) becomes 

[) = Hl{l-'iX + ^l3{l-'i\ + Q^i){l-2,\ + 2/i) Hi 

+97 (1 - 3A + 6/i)^ (5 - 15A + 12^) , (7.9) 

which has the following two non-trivial solutions, 



^° 2(l-3A + 6/u)(5-15A + 12^)7 9 (1 - 3A + 2^)^ /32 )' 

(7.10) 

as long as the r.h.s. is real and positive. If 



(l-3A + 2^)/3 / / 4(1-3A)(5-15A + 12^)7 



4(1-3A) (5-15A + 12^)7 



9 (1 - 3A + 2p) 2/32 



« 1 , (7.11) 



one of the two solutions is much smaller than the other solution. Then one may regard 
that the larger solution corresponds to the inflation in the early universe and the smaller 
one to the late-time acceleration. 

More examples of F (^Rj theory which can contain more than one dS solution, such 
that inflation and dark energy epochs can be explained under the same mechanism Ref. p6| 
may be considered. First of all, as generalization of the model ([7.8|), a general polynomial 
function may be discussed 



m 



F(i?) = ^a„i?", (7.12) 



n=l 



Here an are coupling constants. Using the equation (7^), it yields the algebraic equation. 



n=l n=l 

By a qualitative analysis, one can see that the number of positive real roots, i.e., of the 
de Sitter points, depends completely on the sign of the coupling constants a„. Then, by 
a proper choice, F (^Rj gravity can well explain dark energy and inflationary epochs in 
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a unified natural way. Even it could predict the existence of more than two accelerated 
epochs, which could resolve the coincidence problem. 
Let us now consider an explicit example 

F (R) = - ^ ^ , (7.14) 



where a, /?, 7, n are constants. By introducing this function in (7/7), it is straightforward to 
show that for the function ( [7.14| ), there are several de Sitter solutions. In order to simplify 
this example, let us consider the case n = 2, where the equation ([7. 71 ) yields, 

7 - 37A - 3/3i7^(l - 3A + + 27aH^{-l + 3A - 4/x)(l - 3A + e/x)^ = . (7.15) 

The solutions are given by 

= 7 ^ ^TTT ^ r - 3A + Gfif 

° 18a(l -3A + 6/x)2(-l + 3A-4^) 

±^/{l - 3A + 6/i)2 [12a7(-l + 3A)(-1 + 3A - + /32(1 - 3A + 6;u)2]| . (7.16) 

Then, by a proper choice of the free parameters of the model, two positive roots of the 
equation ( 7.15| ) are solutions. Hence, such a model can explain inflationary and dark energy 
epochs in unified manner. 



8. Discussion 



In summary, in this work we aimed to resolve (at least, partially) the inconsistency problems 
of the projectable HL gravity. First of all, it is demonstrated that some versions of F{R) 
HL gravity may have stable de Sitter solution and unstable flat space solution. As a result, 
the spectrum analysis showing the presence of scalar graviton is not applied. The whole 
spectrum analysis should be redone for de Sitter background. 

Second, U{1) extension of F{R) HL gravity is formulated in two alternative approaches. 
Hamiltonian structure of U{1) invariant F{R) gravity is investigated in all detail. The whole 
constraints system is derived and different particular cases corresponding to conditions for 
derivatives of function F are studied. It is demonstrated that in some cases the Hamiltonian 
structure of the theory coincides with the one of U{1) invariant HL gravity that conforms 
consistency of our approach. The analysis of fluctuations of U{1) invariant F{R) HL gravity 
is performed. It is shown that like in case of U{1) HL gravity the scalar graviton ghost 
does not emerge. This opens good perspectives for consistency of such class of models. It 
is also interesting that spatially-flat FRW equations for U{1) invariant F{R) gravity turn 
out to be just the same as for the one without U{1) symmetry. This indicates that all 
(spatially-flat FRW) cosmological predictions of viable conventional F{R) gravity are just 
the same as for its HL counterpart (with special parameters choice). 
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